**** PROGRAM FILE FOR SI4: Model Sepcification Comparisons ****
 
 
  

use data_final.dta, clear 




**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************



*** DECLARE DATABASE AS A PANEL (AGENCY ID * FISCAL YEAR): 

xtset agenid fyear, yearly

*
*
*


**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************



*** I. ESTIMATE EFFECTS OF BAA REFORMS ONBUDGETARY CONTROL [USING BUDGETARY GROWTH AS DEPENDENT VARIABLE] ***


*** ESTIMATE EXECUTIVE BUDGET CONTROL USING ABADIE'S SEMI-PARAMETRIC PROPENSITY SCORE WEIGHTED ESTIMATOR OF HETEROGENEOUS AVERAGE TREATMENT EFFECTS ***

absdid execbudchangetot, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateect
*
absdid execbudchangetot, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ecti

lincomest _cons + executiveagency
estimate store ecte

absdid execbudchangetot , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateect1

absdid execbudchangetot , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ecti1

lincomest _cons + executiveagency
estimate store ecte1

absdid execbudchangetot , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateect2

absdid execbudchangetot , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ecti2

lincomest _cons + executiveagency
estimate store ecte2

absdid execbudchangetot , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability) sle
estimate store ateect3

absdid execbudchangetot , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability) yxvar(executiveagency) sle
estimate store ecti3

lincomest _cons + executiveagency
estimate store ecte3
*
*
*
absdid execbudchangereg, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateecr
*
absdid execbudchangereg, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ecri

lincomest _cons + executiveagency
estimate store ecre

absdid execbudchangereg , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateecr1

absdid execbudchangereg , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ecri1

lincomest _cons + executiveagency
estimate store ecre1

absdid execbudchangereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateecr2

absdid execbudchangereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ecri2

lincomest _cons + executiveagency
estimate store ecre2

absdid execbudchangereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability) sle
estimate store ateecr3

absdid execbudchangereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability) yxvar(executiveagency) sle
estimate store ecri3

lincomest _cons + executiveagency
estimate store ecre3
*
*
*
*
*** ESTIMATE CONGRESSIONAL BUDGET CONTROL USING ABADIE'S SEMI-PARAMETRIC PROPENSITY SCORE WEIGHTED ESTIMATOR OF HETEROGENEOUS AVERAGE TREATMENT EFFECTS ***

absdid congbudchangetot, tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecct
*
absdid congbudchangetot, tvar(bureauofbudget)  xvar(housemajorityparty  presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ccti

lincomest _cons + executiveagency
estimate store ccte

absdid congbudchangetot , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecct1

absdid congbudchangetot , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ccti1

lincomest _cons + executiveagency
estimate store ccte1

absdid congbudchangetot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecct2

absdid congbudchangetot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ccti2

lincomest _cons + executiveagency
estimate store ccte2

absdid congbudchangetot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest) sle
estimate store atecct3

absdid congbudchangetot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest) yxvar(executiveagency) sle
estimate store ccti3

lincomest _cons + executiveagency
estimate store ccte3
*
*
*
absdid congbudchangereg, tvar(bureauofbudget)  xvar(housemajorityparty  presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateccr

*
absdid congbudchangereg, tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ccri

lincomest _cons + executiveagency
estimate store ccre

absdid congbudchangereg , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateccr1


absdid congbudchangereg , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ccri1

lincomest _cons + executiveagency
estimate store ccre1

absdid congbudchangereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateccr2

absdid congbudchangereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store ccri2

lincomest _cons + executiveagency
estimate store ccre2

absdid congbudchangereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest) sle
estimate store ateccr3

absdid congbudchangereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest) yxvar(executiveagency) sle
estimate store ccri3

lincomest _cons + executiveagency
estimate store ccre3


*** Sensitivity Plots for Budget Control Models***
*** Figure SI-4A. Executive Budget Control: Comparison of Alternative Model Specifications ***

coefplot (ateect, rename(_cons="ATE (Reported Model)")\ ///
			ateect1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			ateect2, rename(_cons="ATE (Alternative Model 2)")\ ///
		 	ateect3, rename(_cons="ATE (Alternative Model 3)")) ///
         (ecte, rename((1)="ATE: E (Reported Model)")\       ///
			ecte1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			ecte2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			ecte3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (ecti, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			ecti1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			ecti2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			ecti3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Total Budget) ///
		  ||                       ///
		 (ateecr, rename(_cons="ATE (Reported Model)")\ ///
			ateecr1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			ateecr2, rename(_cons="ATE (Alternative Model 2)")\ ///
			ateecr3, rename(_cons="ATE (Alternative Model 3)")) ///
         (ecre, rename((1)="ATE: E (Reported Model)")\       ///
			ecre1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			ecre2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			ecre3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (ecri, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			ecri1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			ecri2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			ecri3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Regular Budget) ///
		  ||,                        ///
          nokey xline (0, lcolor(black) lwidth(thin) lpattern(dash)) grid(w) ciopts(recast(rcap)) ///
		  nooffsets msize(medsmall) ///
		  order("ATE (Reported Model)" "ATE (Alternative Model 1)" "ATE (Alternative Model 2)" "ATE (Alternative Model 3)" . ///
                "ATE: E (Reported Model)" "ATE: E (Alternative Model 1)" "ATE: E (Alternative Model 2)" "ATE: E (Alternative Model 3)" . ///
                "ATE: I (Reported Model)" "ATE: I (Alternative Model 1)" "ATE: I (Alternative Model 2)" "ATE: I (Alternative Model 3)" ) ///
		  byopts(note("Point Estimates and Corresponding 95% Confidence Intervals", j(right) place(seast) size(vsmall)) ///
		  title ("Figure SI-4A. 1921 Budget and Accounting Act Average Treatment Effects on " ///
		  "Executive Budget Control: Comparison of Alternative Model Specifications", size (med)span))
		  
graph save "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\gph\Figure SI-4A.gph", replace
graph export "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\Figure SI-4A.png", replace

 
*** Figure SI-4B. Legislative Budget Control: Comparison of Alternative Model Specifications ***

coefplot (atecct, rename(_cons="ATE (Reported Model)")\ ///
			atecct1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			atecct2, rename(_cons="ATE (Alternative Model 2)")\ ///
			atecct3, rename(_cons="ATE (Alternative Model 3)")) ///
         (ccte, rename((1)="ATE: E (Reported Model)")\       ///
			ccte1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			ccte2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			ccte3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (ccti, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			ccti1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			ccti2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			ccti3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Total Budget) ///
		  ||                       ///
		 (ateccr, rename(_cons="ATE (Reported Model)")\ ///
			ateccr1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			ateccr2, rename(_cons="ATE (Alternative Model 2)")\ ///
			ateccr3, rename(_cons="ATE (Alternative Model 3)")) ///
         (ccre, rename((1)="ATE: E (Reported Model)")\       ///
			ccre1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			ccre2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			ccre3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (ccri, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			ccri1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			ccri2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			ccri3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Regular Budget) ///
		  ||,                        ///
          nokey xline (0, lcolor(black) lwidth(thin) lpattern(dash)) grid(w) ciopts(recast(rcap)) ///
		  nooffsets msize(medsmall) ///
		  order("ATE (Reported Model)" "ATE (Alternative Model 1)" "ATE (Alternative Model 2)" "ATE (Alternative Model 3)" . ///
                "ATE: E (Reported Model)" "ATE: E (Alternative Model 1)" "ATE: E (Alternative Model 2)" "ATE: E (Alternative Model 3)" . ///
                "ATE: I (Reported Model)" "ATE: I (Alternative Model 1)" "ATE: I (Alternative Model 2)" "ATE: I (Alternative Model 3)" ) ///
		  byopts(note("Point Estimates and Corresponding 95% Confidence Intervals", j(right) place(seast) size(vsmall)) ///
		  title ("Figure SI-4B. 1921 Budget and Accounting Act Average Treatment Effects on " ///
		  "Legislative Budget Control: Comparison of Alternative Model Specifications", size (med)span))
		  
graph save "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\gph\Figure SI-4B.gph", replace
graph export "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\Figure SI-4B.png", replace


**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************




*** II. ESTIMATE EFFECTS OF BAA REFORMS ON BUDGETARY STABILITY IN NATURAL LOGARITHMS AS DEPENDENT VARIABLE ***



*** ESTIMATE EXECUTIVE BUDGET STABILITY USING ABADIE'S SEMI-PARAMETRIC PROPENSITY SCORE WEIGHTED ESTIMATOR OF HETEROGENEOUS AVERAGE TREATMENT EFFECTS ***

absdid lnexecbudstabletot, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateest
*
absdid lnexecbudstabletot, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store esti

lincomest _cons + executiveagency
estimate store este

absdid lnexecbudstabletot , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateest1

absdid lnexecbudstabletot , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store esti1

lincomest _cons + executiveagency
estimate store este1

absdid lnexecbudstabletot , tvar(bureauofbudget)  xvar(presidentparty presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateest2

absdid lnexecbudstabletot , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store esti2

lincomest _cons + executiveagency
estimate store este2

absdid lnexecbudstabletot , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability) sle
estimate store ateest3

absdid lnexecbudstabletot , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability)  yxvar(executiveagency) sle
estimate store esti3

lincomest _cons + executiveagency
estimate store este3
*
*
*
absdid lnexecbudstablereg, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateesr
*
absdid lnexecbudstablereg, tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store esri

lincomest _cons + executiveagency
estimate store esre

absdid lnexecbudstablereg , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateesr1

absdid lnexecbudstablereg , tvar(bureauofbudget)  xvar(agencyheadstability realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store esri1

lincomest _cons + executiveagency
estimate store esre1

absdid lnexecbudstablereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store ateesr2

absdid lnexecbudstablereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store esri2

lincomest _cons + executiveagency
estimate store esre2

absdid lnexecbudstablereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability) sle
estimate store ateesr3

absdid lnexecbudstablereg , tvar(bureauofbudget)  xvar(presidentparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
agencyheadstability)  yxvar(executiveagency) sle
estimate store esri3

lincomest _cons + executiveagency
estimate store esre3


*
*
*

*** ESTIMATE CONGRESSIONAL BUDGET STABILITY USING ABADIE'S SEMI-PARAMETRIC PROPENSITY SCORE WEIGHTED ESTIMATOR OF HETEROGENEOUS AVERAGE TREATMENT EFFECTS ***

absdid lncongbudstabletot, tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecst
*
absdid lncongbudstabletot, tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store csti

lincomest _cons + executiveagency
estimate store cste

absdid lncongbudstabletot , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecst1

absdid lncongbudstabletot , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store csti1

lincomest _cons + executiveagency
estimate store cste1

absdid lncongbudstabletot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecst2

absdid lncongbudstabletot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store csti2

lincomest _cons + executiveagency
estimate store cste2

absdid lncongbudstabletot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest) sle
estimate store atecst3

absdid lncongbudstabletot , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest)  yxvar(executiveagency) sle
estimate store csti3

lincomest _cons + executiveagency
estimate store cste3
*
*
*
absdid lncongbudstablereg, tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecsr
*
absdid lncongbudstablereg, tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store csri

lincomest _cons + executiveagency
estimate store csre

absdid lncongbudstablereg , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecsr1

absdid lncongbudstablereg , tvar(bureauofbudget)  xvar(appropcommchairstabilityest realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store csri1

lincomest _cons + executiveagency
estimate store csre1

absdid lncongbudstablereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atecsr2

absdid lncongbudstablereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
realgnpgrowth realfeddeficitgnp lnagencyage)  yxvar(executiveagency) sle
estimate store csri2

lincomest _cons + executiveagency
estimate store csre2

absdid lncongbudstablereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest) sle
estimate store atecsr3

absdid lncongbudstablereg , tvar(bureauofbudget)  xvar(housemajorityparty presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
appropcommchairstabilityest)  yxvar(executiveagency) sle
estimate store csri3

lincomest _cons + executiveagency
estimate store csre3

*** Sensitivity Plots for Budget Stability Models***
*** Figure SI-4C. Executive Budget Stability: Comparison of Alternative Model Specifications ***

coefplot (ateest, rename(_cons="ATE (Reported Model)")\ ///
			ateest1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			ateest2, rename(_cons="ATE (Alternative Model 2)")\ ///
			ateest3, rename(_cons="ATE (Alternative Model 3)")) ///
         (este, rename((1)="ATE: E (Reported Model)")\       ///
			este1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			este2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			este3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (esti, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			esti1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			esti2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			esti3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Total Budget) ///
		  ||                       ///
		 (ateesr, rename(_cons="ATE (Reported Model)")\ ///
			ateesr1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			ateesr2, rename(_cons="ATE (Alternative Model 2)")\ ///
			ateesr3, rename(_cons="ATE (Alternative Model 3)")) ///
         (esre, rename((1)="ATE: E (Reported Model)")\       ///
			esre1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			esre2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			esre3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (esri, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			esri1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			esri2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			esri3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Regular Budget) ///
		  ||,                        ///
          nokey xline (0, lcolor(black) lwidth(thin) lpattern(dash)) grid(w) ciopts(recast(rcap)) ///
		  nooffsets msize(medsmall) ///
		  order("ATE (Reported Model)" "ATE (Alternative Model 1)" "ATE (Alternative Model 2)" "ATE (Alternative Model 3)" . ///
                "ATE: E (Reported Model)" "ATE: E (Alternative Model 1)" "ATE: E (Alternative Model 2)" "ATE: E (Alternative Model 3)" . ///
                "ATE: I (Reported Model)" "ATE: I (Alternative Model 1)" "ATE: I (Alternative Model 2)" "ATE: I (Alternative Model 3)" ) ///
		  byopts(note("Point Estimates and Corresponding 95% Confidence Intervals", j(right) place(seast) size(vsmall)) ///
		  title ("Figure SI-4C. 1921 Budget and Accounting Act Average Treatment Effects on " ///
		  "Executive Budget Stability: Comparison of Alternative Model Specifications", size (med)span))
		  
graph save "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\gph\Figure SI-4C.gph", replace
graph export "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\Figure SI-4C.png", replace

*** Figure SI-4D. Legislative Budget Stability: Comparison of Alternative Model Specifications ***

coefplot (atecst, rename(_cons="ATE (Reported Model)")\ ///
			atecst1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			atecst2, rename(_cons="ATE (Alternative Model 2)")\ ///
			atecst3, rename(_cons="ATE (Alternative Model 3)")) ///
         (cste, rename((1)="ATE: E (Reported Model)")\       ///
			cste1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			cste2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			cste3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (csti, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			csti1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			csti2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			csti3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Total Budget) ///
		  ||                       ///
		 (atecsr, rename(_cons="ATE (Reported Model)")\ ///
			atecsr1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			atecsr2, rename(_cons="ATE (Alternative Model 2)")\ ///
			atecsr3, rename(_cons="ATE (Alternative Model 3)")) ///
         (csre, rename((1)="ATE: E (Reported Model)")\       ///
			csre1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			csre2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			csre3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (csri, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			csri1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			csri2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			csri3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Regular Budget) ///
		  ||,                        ///
          nokey xline (0, lcolor(black) lwidth(thin) lpattern(dash)) grid(w) ciopts(recast(rcap)) ///
		  nooffsets msize(medsmall) ///
		  order("ATE (Reported Model)" "ATE (Alternative Model 1)" "ATE (Alternative Model 2)" "ATE (Alternative Model 3)" . ///
                "ATE: E (Reported Model)" "ATE: E (Alternative Model 1)" "ATE: E (Alternative Model 2)" "ATE: E (Alternative Model 3)" . ///
                "ATE: I (Reported Model)" "ATE: I (Alternative Model 1)" "ATE: I (Alternative Model 2)" "ATE: I (Alternative Model 3)" ) ///
		  byopts(note("Point Estimates and Corresponding 95% Confidence Intervals", j(right) place(seast) size(vsmall)) ///
		  title ("Figure SI-4D. 1921 Budget and Accounting Act Average Treatment Effects on " ///
		  "Legislative Budget Stability: Comparison of Alternative Model Specifications", size (med)span))
		  
graph save "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\gph\Figure SI-4D.gph", replace
graph export "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\Figure SI-4D.png", replace



**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************



*** III. ESTIMATE EFFECTS OF BAA REFORMS ON EXECUTIVE-LEGISLATIVE BDUGETARY COHERENCE IN NATURAL LOGARITHMS (+1) AS DEPENDENT VARIABLE ***



*** ESTIMATE EXECUTIVE-CONGRESSIONAL BUDGET COHERENCE USING ABADIE'S SEMI-PARAMETRIC PROPENSITY SCORE WEIGHTED ESTIMATOR OF HETEROGENEOUS AVERAGE TREATMENT EFFECTS ***

absdid ln1execbudcoherencetot, tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) sle  
estimate store atebct
*
absdid ln1execbudcoherencetot, tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store bcti

lincomest _cons + executiveagency
estimate store bcte

absdid ln1execbudcoherencetot , tvar(bureauofbudget)  xvar(sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) sle 
estimate store atebct1

absdid ln1execbudcoherencetot , tvar(bureauofbudget)  xvar(sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store bcti1

lincomest _cons + executiveagency
estimate store bcte1

absdid ln1execbudcoherencetot , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atebct2

absdid ln1execbudcoherencetot , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store bcti2

lincomest _cons + executiveagency
estimate store bcte2

absdid ln1execbudcoherencetot , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest) sle 
estimate store atebct3

absdid ln1execbudcoherencetot , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest) yxvar(executiveagency) sle
estimate store bcti3

lincomest _cons + executiveagency
estimate store bcte3
*
*
*
absdid ln1execbudcoherencereg, tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atebcr
*
absdid ln1execbudcoherencereg, tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store bcri

lincomest _cons + executiveagency
estimate store bcre

absdid ln1execbudcoherencereg , tvar(bureauofbudget)  xvar(sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atebcr1

absdid ln1execbudcoherencereg , tvar(bureauofbudget)  xvar(sharedtenureest realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store bcri1

lincomest _cons + executiveagency
estimate store bcre1

absdid ln1execbudcoherencereg , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) sle
estimate store atebcr2

absdid ln1execbudcoherencereg , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt realgnpgrowth realfeddeficitgnp lnagencyage) yxvar(executiveagency) sle
estimate store bcri2

lincomest _cons + executiveagency
estimate store bcre2

absdid ln1execbudcoherencereg , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest) sle
estimate store atebcr3

absdid ln1execbudcoherencereg , tvar(bureauofbudget)  xvar(presidentialstabilityadministrat housemajoritypartystability senatemajoritypartystability ///
unifiedpartygovt sharedtenureest) yxvar(executiveagency) sle
estimate store bcri3

lincomest _cons + executiveagency
estimate store bcre3
*
*
*
*** Sensitivity Plots for Budgetary Coherence Models***
*** Figure SI-4E. Budget Coherence: Comparison of Alternative Model Specifications ***

coefplot (atebct, rename(_cons="ATE (Reported Model)")\ ///
			atebct1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			atebct2, rename(_cons="ATE (Alternative Model 2)")\ ///
			atebct3, rename(_cons="ATE (Alternative Model 3)")) ///
         (bcte, rename((1)="ATE: E (Reported Model)")\       ///
			bcte1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			bcte2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			bcte3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (bcti, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			bcti1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			bcti2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			bcti3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Total Budget) ///
		  ||                       ///
		 (atebcr, rename(_cons="ATE (Reported Model)")\ ///
			atebcr1, rename(_cons="ATE (Alternative Model 1)")\ /// 
			atebcr2, rename(_cons="ATE (Alternative Model 2)")\ ///
			atebcr3, rename(_cons="ATE (Alternative Model 3)")) ///
         (bcre, rename((1)="ATE: E (Reported Model)")\       ///
			bcre1, rename((1)="ATE: E (Alternative Model 1)")\ ///
			bcre2, rename((1)="ATE: E (Alternative Model 2)")\ ///
			bcre3, rename((1)="ATE: E (Alternative Model 3)")) ///
         (bcri, drop(executiveagency) rename(_cons="ATE: I (Reported Model)")\ ///
			bcri1, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 1)")\ ///
			bcri2, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 2)")\ ///
			bcri3, drop(executiveagency) rename(_cons="ATE: I (Alternative Model 3)")) ///
          ,  bylabel(Regular Budget) ///
		  ||,                        ///
          nokey xline (0, lcolor(black) lwidth(thin) lpattern(dash)) grid(w) ciopts(recast(rcap)) ///
		  nooffsets msize(medsmall) ///
		  order("ATE (Reported Model)" "ATE (Alternative Model 1)" "ATE (Alternative Model 2)" "ATE (Alternative Model 3)" . ///
                "ATE: E (Reported Model)" "ATE: E (Alternative Model 1)" "ATE: E (Alternative Model 2)" "ATE: E (Alternative Model 3)" . ///
                "ATE: I (Reported Model)" "ATE: I (Alternative Model 1)" "ATE: I (Alternative Model 2)" "ATE: I (Alternative Model 3)" ) ///
		  byopts(note("Point Estimates and Corresponding 95% Confidence Intervals", j(right) place(seast) size(vsmall)) ///
		  title ("Figure SI-4E. 1921 Budget and Accounting Act Average Treatment Effects on " ///
		  "Budget Coherence: Comparison of Alternative Model Specifications", size (med)span))
		  
graph save "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\gph\Figure SI-4E.gph", replace
graph export "C:\Users\Jin\Dropbox\BACA paper\Graphs\Governance\Figure SI-4E.png", replace


*
*
*


**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************
**********************************************************************************************************************************************************************************************************




